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A Dusty Component to the Gaseous Debris Disk around the 
White Dwarf SDSS J1228+1040 

C. S. Brinkworth/ B. T. Gansicke,^ T. R. Marsh,^ D. W. Hoard, ^ C. Tappert,^ 

ABSTRACT 

We present ISAAC spectroscopy and ISAAC, UKIDSS and Spitzer Space Tele- 
scope broad-band photometry of SDSS J1228-I-1040 - a white dwarf for which 
evidence of a gaseous metal-rich circumstellar disk has previously been found 
from optical emission lines. The data show a clear excess in the near- and mid- 
infrared, providing compelling evidence for the presence of dust in addition to 
the previously identified gaseous debris disk around the star. The infrared excess 
can be modelled in terms of an optically thick but geometrically thin disk. We 
find that the inner disk temperatures must be relatively high (~1700 K) in order 
to fit the SED in the near-infrared. These data provide the first evidence for the 
co-existence of both gas and dust in a disk around a white dwarf, and show that 
their presence is possible even around moderately hot (~22,000 K) stars. 

Subject headings: Stars: individual (SDSS J1228+1040) — circumstellar matter 
- white dwarfs - infrared: stars 



1. Introduction 

Despite the numerous extrasolar planets discovered over the last decade, the fate of 
these systems during the final stages of the parent star's evolution is still uncertain. A 
clue to their destiny was discovered during a search for cool companions to white dwarfs . 



G29-38 was found to display a strong near- infrared excess (jZuckerman fc BecklinI 119871 ). 
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but deep imaging and asteroseismological studies of the system ruled out the presence of 
a brown dwarf companion. Mid-infrared photometry showed that this excess is even more 
pronounced at 10 fim suggesting that the emission is due to du st around the white dwarf 
( IGraham et aLl Il990l : iTelescd iJoy fc SiskI : iTokunaga et aLlllQQOl ). Spitzer Space Telescope 
observations spectr oscopically confirmed the presence of circumstellar dust around G29-38 



von Hippel et ah 



( Reach et al.l 12005), wh ich has since been modelled as a flat, opaque, circumstellar disk 



20071 ). The atmosphere of G29-38 



las been found to be enriched with 



metals, implying that the star is currently accreting (IKoester et al.l 119971 ). probably from 
this dusty disk. 

The more likely origin for the disk is the tidal disruption of either comets or asteroids 
( iDebes fc SigurdssonI l2002l : I Jural l2003l ). with asteroids the more probable candidates, as 
they can explain the large amount of metals accreted by the white dwarfs. Comets are 
also dynamically less favourable, g i ven th eir eccentric orbits and long orbital periods. The 
approach used by lvon Hippel et al.l (120071 ) to model the disk was ba sed on a simple, optically 
thick, geometrically thin dusty disk model described in I Jural (120031 ). 



A further 11 white dwarfs have been found to disp l ay a similar infrare d excess (iBecklin et al 



2005 



2008a 



Kihc et all |2005|, l2006l : Ivon Hippel et all l2007l : lJura et al.l l2007bl . l2008l : iFarihi et al. 



bj), most of which are cool, metal- rich DAZ white dwarfs, with effective temperatures 
Tcff < 15000 K (with the exceptio ns of PG 1()15+16 1, which has Teg ~19300 K, and PG 1457- 
086 with Teff ~ 20400 K), leading iKilic et al.l (120061 ) to suggest that du sty disks around white 



dwar f s shou ld be sublimated for Tgfj > 16000 - 20000 K. More recently, iGansicke et al.l (12006 



20071 . 12OO8I ) identified three more moderately hot white dwarfs showing double-peaked cal- 
cium lines indicative of a rotating circumstellar debris disk - SDSS J122859. 93-1-104032. 9 
(hereafter SDSS 1228+1040), SDSS J104341. 53+085558.2 and SDSS J084539. 17+225728.0. 

Here we report the detection of a significant i nfrared excess in one of these systems. 



SDSS 1228+1040. The white dwarf was identified bvlGansicke et al.l ( 20061) in Data Release 



6 of the Sloan Digital Sky Survey (lAdelman-McCarthy et al 



200 



i 



Gansicke et all (l2006f ) 



noted very unusual emission lines of the Ca II 850-866 nm triplet, as well as weaker emis- 
sion lines of Fe II at 502 nm and 517 nm. The line profiles of the Ca II triplet display a 
double-peaked morphology, indicative of a gaseous, rotating disk structure. Time-resolved 
spectroscopy ruled out any significant radial velocity variations, and the absence of Balmer 
and helium emission lines implied that the gaseous disk around the white dwarf must be 
extremely deficient in volatile materials, ruling out the presence of a stellar binary com- 
panion. The velocity of the Ca II line peaks indicate that the outer radius of the gaseous 
circumstellar disk is ~ 1.2 R©. 
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2. Observations and Data Reduction 

2.1. Spitzer Space Telescope data 

We were awarded 2.4 hours of Spitzer Space Telescope time in Cycle-4 (P40048), plus 
3.0 hours in Cycle-5 (P50118) to search for a dusty component to the gaseous circumstel- 
lar disk around SPSS 1228+10 40 (see Table[I] for a full list of observations). We obtained 



band Imaging Photometer (MIPS; iRieke et al. 



Infrared Array Camera (IR AC: iFazio et 3.1.1120041) data in all four channels (3.6 - 8.0 yiim), 



Infrared Spectrograph (IRS: lHouck et al.ll2004 ) Blue Peak-Up imaging at 16yum, and Multi- 



20041 ) imaging at 24/im. 



The IRAC data reduction was carried out on the Basic Calibrated Data (BCD) frames as 
follows. We corrected all of the BCDs for array-location-dependence using the corresponding 
correction frames provided by the Spitzer Science Center (SSC). These were applied after 
uncorrecting the array-location-dependence frames for pixel distortion in order to prevent a 
repeat correction for pixel distortion when creating the mosaic. The array-loc ation-corrected 



BCDs were then all combined using the standard SSC software, MOPEX (IMakovoz et al. 



20061 . MOsaicker and Point source Extraction) with dual-outlier rejection to create a single 
mosaicked image for each channel. The mosaic pixel scales were set to the default 1.22"/pix. 
As a check on our reduction method, we also split the BCDs into batches of 7 frames, creating 
a mosaic from each batch in order to enable us to test the repeatability of the photometry 
and calculate an independent estimate of our uncertainties. 

The IRS Peak-Up imaging frames were combined in much the same way, using MOPEX 
to combine the BCDs with dual-outlier rejection. Again, we created a mosaic from both the 
entire data set, and from sub-sets of the data, so that we could test the repeatability of the 
photometry and derive realistic photometric uncertainties from our data. The pixel scales 
for the mosaics were set to the default 1.8"/pix. A similar method was applied to the MIPS 
data, where we created a separate mosaic for each of the Exposure ID sets (35 BCDs in each 
set). This gave us 21 mosaics in total, all of which were created using multiframe temporal 
outlier rejection and the default pixel scale of 2.45"/pix. The mosaicked images for all of the 
Spitzer data, with SDSS 1228-1-1040 marked, can be found in FigurelH 

We carried out aperture photometry on our IRAC and IRS PU mosaics using IRAF, 
using apertures of 3.0 pixels for the IRAC mosaics and 2.0 pixels for the IRS mosaics. 
The sky subtraction was applied using an annulus of 10-20 pixels in radius. In IRAC 
channels 1 and 2 there were background sources included in the sky annuli. IRAF rejects 
bright pixels in the sky annuli, but we carried out a double-check of the sky levels using 
nearby apertures, and found them to be virtually identical to the levels originally measured. 
The photometry was then converted from MJy/sr to /iJy and aperture-corrected to the 
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calibration aperture sizes using the standard aperture corrections provided by the SSC. The 
IRAC channel 1 photometry was further corrected for pixel-phase dependence in accordance 
with SSC recommendations. No color correction was applied since we quote the isophotal 
effective wavelengths, thus eliminating most of the color-dependency of the flux calibration 
and rendering the magnitude of the correction negligible compared to our uncertainties. 

We derived our uncertainty estimates for the IRAC and IRS data using the stan- 
dard deviation of the scatter in the sub-set mosaics at each wavelength. These were then 
added in quadratu re to the absolute calibration estimates of the IRAC and IRS instruments 
faeach et al.lboosl . SSC websitfl). 



Looking at the MIPS mosaics, we found that the target is very faint, and there is bright 
source within ~ 5 — 6 pixels. While it is tempting to suspect that the bright source may, 
in fact, be the target, the pointing of Spitzer is good to within 0.5", while the pixel scale 
of the mosaic is 2.45"/pix, and so there is no doubt that the fainter of the two sources 
is SDSS 1228+1040. This is corroborated by matching the positions with respect to the 
surrounding stars. We used MOPEX to carry out a source extraction and subtraction on 
each mosaic, removing all of the high SNR sources from the field (including the close, bright 
neighbour), before carrying out IRAF aperture photometry on our lower SNR, unsubtracted 
target. We used an aperture of 3"(1.224 pixels) to extract the photometry, with a background 
annulus of 20-32", and applied the standard aperture corrections as given in the MIPS Data 
Handbook. We double-checked the flux density we obtained by performing MOPEX PRF- 
fitting on our target in the original mosaics. The two methods agreed to within 10%, while the 
standard deviation in the spread in flux densities from the 21 mosaics gave us uncertainties 
of ~50%. The final flux density was taken as the unweighted median in the 21 measurements 
from the aperture photometry. 



2.2. ISAAC data 

The infrared spectroscopic data were taken with the Infrared Spectrometer And Array 
Camera (ISAAC) mounted at Antu (UTl), VLT, Paranal, Chile. The short-wavelength 
grism was used in its low-resolution mode at central wavelengths 1.06 /im (z), 1.25 /im (J), 
1.65 /im (H), and 2.20 /xm (K). An 0.8" slit width resulted in resolving powers in the range 
650-750. The observations were conducted in service mode on two nights, 2007 April 03 (J 
and H) and 2007 April 07 {z and K). Flat fields and wavelength calibration with a Xe-Ar 
lamp (Ar only for the z spectra) were taken at the start of the night, and telluric standards 



http:/ /ssc. spitzer. caltech.edu/irac/documents/reach2005.pdf 
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were observed within 1 h and an airmass difference AM{z) = 0.2 of tfie target spectra. The 
data were obtained as sets of 22 spectra per central wavelength with individual exposure 
times of 100 s. The individual spectra were taken at different positions on the detector (AB 
cycles plus jittering) to minimize the effect of bad pixels and cosmic rays. 

For the reduction, IRAF packages were applied. As first steps, the data were flatfielded 
and "straightened" via a two-dimensional wavelength calibration, thus correcting for the 
positional dependence of the dispersion. Corresponding AB pairs were subtracted from each 
other, and the resulting data then were combined to a single image. Subsequently, the target 
and standard spectra were extracted. In the case of the z spectrum it proved ne cessary to 



adjust the wavelength calibration using the night sky lines (IRousselot et al.ll2000l ). 



The following stars were used as telluric standards: HD 106807 (AlV, T^s = 9 400K) 
for the z data. Hip 43699 (B4V, 17 000 K) for J and H, and Hip 88857 (B3V, 19 000 K) for 
K. The intrinsic lines of the standards were fitted with a Voigt profile and subtracted from 
the infrared spectrum, to yield a pure atmospheric absorption spectrum. For the z data, the 
intrinsic lines were submerged in atmospheric absorption bands and could not be corrected 
for. The target spectra were then divided by their corresponding telluric spectra, after the 
latter had been shifted and scaled to match the position and depths of the atmospheric fea- 
tures in the target data. Finally, the spectra were multiplied with a blackbody corresponding 
to the effective temperature of the telluric standard, thus recovering the intrinsic SED. 

We used the ISAAC JHK acquisition images to perform differential photometry against 
nearby bright stars that have reliable magnitudes in 2MASS, where we assumed a conserva- 
tive 10 % error for the IR fluxes derived from these measurements. The JHK magnitudes of 
SDSS 1228-1-1040 obtained in this way were then used to correct the flux level of the ISAAC 
spectra for slit losses. 



2.3. UKIDSS data 



SDSS 1228-1-1040 is within the footprint of the Large Area Survey within the third Data 
Release (DR3) of the UKIRT Infrared Deep Sky Survey Jhawrence et al.ll2007l. UK IDSS) . 
We obtained Y JHK aperture magnitudes from the WFCAM archive (IHambly et al.ll2008l ). 
and converted them to fluxes at the effective wavelengths using the zero-point definitions in 
Hewett et al.l (120061 ). The agreement between the UKIDSS and ISAAC JHK fluxes is very 
good. 
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3. Results and Modelling 



The results of the Spitzer, UKIDSS and ISAAC photometry are given in Table[2l along 
with their uncertainties. Figure[2] shows the Spitzer and ISAAC photometry, and the ISAAC 



spectr oscopy plotted on the same axes as the UV and optical data from iGansicke et al. 



(120061 ) ■ The model of the expected flux from the white dwarf (taken from the same paper) 
is also shown. There is an substantial flux excess over the expected emission from the white 
dwarf, starting in the K band and increasing out to 16/im, before dropping at 24/im. 

We attempted to model this excess with a simple optically thick, geometrically thin 
dust di sk model w it h a temperature pro file T nisk oc r~^, with (3 = 3/ 4, as previously used 



e.g. by I Jural (120031 ): iBecklin et al.l (120051 ) and iBrinkworth et al.l (120071 ). In principle, P can 
take values different from 3/4 depending on the size of the dust grains, as the grains will 
cool inefficently once that the grain size is comparable or smaller than the peak wavelength 
of the Plan ck function at the grairi temperature, rgrain ^ b/T (with b the Wien constant, for 
details, see IBrinkworth et al.ll2007l ). In practice, howe ver, we found that value s significantly 
different from 3/4 do not produce viable results. As in IBrinkworth et al.l (120071 ). we therefore 
adopt this value for our modelling. The input parameters are given in Table[3l 

The model is shown in Figure[31 While it fits the near-IR, IRAC and MIPS points well, 
the 16 micron flux density is 2.0 sigma greater than predicted. This 16 /im excess may be 
explained by the presence of a sigi iificant silicate ern ision feature at lO/^m, as seen in the 
similar dusty white dwarf CD 362. lJura et al.l (l2007al ) found that this emission feature has 
a strong red wing that, if present in SDSS 1228+1040, would contribute significantly to the 
IRS Peak-Up band (which has a passband from 13 - 18.5 /im). It is possible that such a 
feature would also affect the 8 /im flux density, but the blue wing of the silicate feature in 
GD362 was less pronounced than the red wing, so the effect is expected to be minor (of 
order a few %). This is not currently accounted for in the model, but could explain the 
slight rise seen at 8 /im from the rest of the IRAC bands. 



4. Discussion and Conclusions 



The ISAAC, UKIDSS and Spitzer Space Telescope data show a noticeable infrared excess 
in SDSS 1228+1040, in excess of the flux density expected from the white dwarf alone. The 
observed spectral energy distribution provides clear and compelhng evidence for the presence 
of a substanti al dusty compo r ient t o the debris disk, in addition to the gaseous component 
discovered by iGansicke et al.l (120061 ) , indicating that both gaseous and dusty debris material 
can coexist around a white dwarf. The presence of dust around SDSS 1228+1040 implies 
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that debris disks ca n survive a r ound r elatively hot parent wh ite dwarfs, contrary to the 
earlier suggestion by lKilic et al.l (120061 ): Ivon Hippel et al.l (120071 ). based on a smaller sample 
of white dwarfs with debris disks, that white dwarfs hotter than around 15,000 K should 
sublimate any orbiting dust disk. 

Our model shows that an optically thick dust disk model can reproduce the SED seen in 
SDSS 1228+1040, except for the 16 fim point, which we currently attribute to the red wing 
of a strong silicate feature at 10/im, as seen in similar systems (e.g. GD 362). We calculated 
a lower-limit estimate of the r nass of dust in the opt ically thick disk using the generic mass 
absorption coefficient used by lBeckwith et al.l (Il990l ) of: 



where we follow their example and adopt = 1 (see iBrinkworth et al.ll2007l . for a more 
in-depth discussion). This leads to a disk mass of at least 1 xlO^^ g, equivalent to ~1.3 
xIO'^^Mmooti or ~0.01 times the mass of Ceres. 

There are, however, two main issues with our disk model. Firstly, in order to match 
the rise from the K band to the IRAC channel 1 point at 3.6/xm, we find that the inner 
temperature of the dust disk has to be around 1670 K, which is hotter than the expected 
sublimation temperature of the dust grains (ar ound 1200-1400 K). M oreover, if we esti- 
mate the inner radius of the gaseous disk from iGansicke et al.l (120061 ) using i?inner disk ~ 
(Vpeak/Kiax)^ X -Router disk, whcrc l^peak IS the velocity of the line peaks and Knax is the maxi- 
mum velocity at which emission is seen, then we find that the inner radius of the dusty disk 
at this temperature is smaller (18 Rwd) than the inner radius of the gaseous disk (~27Rwd)- 
It is possible that the inclusion of a 2- or 3-zone, or warped disk model would allow us to fit 
the data using a cooler inner disk temperature (and therefore larger inner dust disk radius) 
but the lack of data points in the mid-IR would make it impossible to usefully constrain such 
a model, and so we do not pursue it further. 

Secondly, we find that the implied outer radius of the dusty disk (~ 107 -Rwd) coincides 
almost per fectly with the outer radius of the optically thin gaseous disk (~ 108i?wD)5 as 
derived by lGansicke et al.l (120061 ). so the gaseous dust disk permeates throughout the dusty 
component. This is not unexpected in itself, as the sublimating dust will likely replenish the 
gas disk. It is possible, however, that the lifetime of t he dust grains will be relatively short 
due to the increased viscosity. iGansicke et al.l (120071 ) determined the photospheric abun- 
dance of Mg in SDSS122 8+1040 to be log(Mg/H) = -4.58 ± 0.06. Using the prescription of 
Koester fc WilkenI (120061 ). and assuming that hydrogen fraction in the circumstellar material 
is depleted by a factor one hundred with respect to its solar value, we estimate an accretion 
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rate of M ~ 3 x lO^gs ^. This is within the range of accretion rates obtained by lFarihi et al. 



(l2008bl ) following a similar approach for the currently known sample of metal-polluted white 
dwarfs with dusty debris disks. We stress, however, that this estimate of the accretion rate 
should is likely to be a strict upper limit, as M is inverse proportional to the hydrogen 
depletion ir i the circumstellar material, and the factor 100 assumed above is probably under- 
estimated (IFriedrich et al.lll999l : iGansicke et al.l 120081 ) . Combining the accretion rate with 
the estimate of the disk mass from above, we estimate a lower limit on the life time of the 
disk of ~ 1000 yr. 



We thank the participants of the Winter AAS 2008 White Dwarf Splinter Session for 
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Fig. 1. — Mosaicked images for SDSS 1228+1040, North up, East left. Top row, left to right: 
IRAC-1 (3.6/im), IRAC-2 (4.5/im), IRAC-3 (5.8/im). Bottow row, left to right: IRAC- 
4 (8.0/xm), IRS PU (16/im), MIPS-24 (24/im). Pixel scales are 1.22'7pix, 1.8'7pix and 
2.45"/pix for the IRAC, IRS PU and MIPS-24 mosaics, respectively. The telescope pointing 
is good to < 0.5" 
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Fig. 2. — Spectral Energy Distribution (SED) for SDSS 1228+1040 from the ultraviolet 
to the mid-infrared. UV spectra and the white dwarf model (dashed line) are taken from 
Gansicke et al.l (120061 ). The near-IR spectroscopy and YJHK photometry were taken with 
ISAAC/VLT and UKIDSS. The mid-IR observations from 3.6 - 24/im were taken with the 
Spitzer Space Telescope IRAC, IRS Blue Peak-Up and MIPS instruments. The SED shows 
a significant flux density excess from the K-band to longer wavelengths, compared to the 
white dwarf alone. 
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Fig. 3. — The infrared spectral energy distribution of SDSS 1228+1040, compared with a 
model consisting of a white dwarf (dashed line) and an optically thick, geometrically thin 
circumstellar dust disk (dash- dot- dot- dot line). The significant excess at 16/im may be due 
to a strong silicate emission feature. For more details see Section[31 
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Table 1. Observation log for Spizer Space Telescope data 



Instrument/ 


Wavelength 


AOR Key 


Total Integration Time 


Date 


Pipelin 


Channel 


microns 




(s) 






IRAC Ch. 1 


5.550 


22247936 


1000 


2007-06-30 


S16.1 


IRAC Ch. 2 


4.493 


22247936 


1000 


2007-06-30 


S16.1 


IRAC Ch. 3 


5.731 


22247936 


1000 


2007-06-30 


S16.1 


IRAC Ch. 4 


7.872 


22247936 


1000 


2007-06-30 


S16.1 


IRS Blue Peak-Up 


15.6 


22248192 


3450 


2007-06-16 


S16.1 


MIPS-24 


23.68 


25459200 


3012 


2008-07-25 


S18.1 


MIPS-24 


23.68 


25459456 


3012 


2008-07-25 


S18.1 


MIPS-24 


23.68 


25459712 


3012 


2008-07-25 


S18.1 
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Table 2. Photometry for SDSS 1228+1040 from Spitzer and ISAAC 



Instrument / 


Wavelength 


Flux Density 


Uncertainty 


Band 


/im 


mJy 


% 


ISAAC / J 


1.235 


0.277 


10 


ISAAC / H 


1.662 


0.193 


10 


ISAAC / K 


2.159 


0.180 


10 


UKIDSS / y 


1.031 


0.363 


1.1 


UKIDSS / J 


1.248 


0.266 


1.5 


UKIDSS / H 


1.631 


0.183 


2.2 


UKIDSS / K 


2.101 


0.167 


3.6 


IRAC / 1 


3.550 


0.235 


4.5 


IRAC / 2 


4.493 


0.235 


4.1 


IRAC / 3 


5.731 


0.232 


7.5 


IRAC / 4 


7.872 


0.244 


12.8 


IRS PU / 16 


15.6 


0.415 


23.0 


MIPS / 24 


23.68 


0.080 


49 
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Table 3. Parameters for the dust disk model discussed in Section[3l 



Inner Disk Temperature (K) 


1670 


Inner Disk Radius (Rwd) 


18 


Outer Disk Temperature (K) 


450 


Outer Disk Radius (Rwd) 


107 


Disk Inclination (degrees) 


70 



